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Abstract 

We study the possible existence of deconfined quark matter in the 
interior of neutron stars using the Nambu-Jona-Lasinio model to de- 
scribe the quark phase. We find that typical neutron stars with masses 
around 1.4 solar masses do not possess any deconfined quark matter in 
their center. This can be traced back to the property of the NJL model 
which suggests a large constituent strange quark mass over a wide range 
of densities. 
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1 Introduction 



At large temperatures or large densities hadronic matter is expected to undergo 
two phase transitions: one which deconfines quarks (and gluons) and one which 
restores chiral symmetry. Up to now it is an unsettled issue whether these 
two phase transitions are distinct or coincide. The more, it is even unclear 
whether there are real phase transitions or only rapid crossover transitions. 
Such transitions have received much attention in heavy ion physics as well as 
in the context of neutron stars which provide a unique environment to study 
cold matter at supernuclear densities Jl|, 0. Even though a deconfmement 
phase transition seems intuitively evident at large enough densities, from a 
theoretical point of view a confirmation of the existence of a deconfined quark 
phase in neutron stars is so far limited by the uncertainties in modeling QCD 
at large densities. All the more it is important to study and compare different 
available models to shed some light on similarities and differences with respect 
to the behavior of matter at large densities as well as on the corresponding 
predictions of neutron star properties like e.g. its mass and radius. In the future 
such experience may prove to be useful if either an improved understanding of 
matter under extreme conditions provides a more exclusive selection between 
the various models or new experimental results on neutron star properties are 
available to set more stringent constraints. 

Usually the quark matter phase is modeled in the context of the MIT bag 
model 0, [| |U as a Fermi gas of u, d, and s quarks. In this model the phe- 
nomenological bag constant -Bmit is introduced to mimic QCD interactions to 
a certain degree. The investigation of such a phase was furthermore stimulated 
by the idea that a quark matter phase composed of almost an equal amount 
of the three lightest quark flavors could be the ground state of nuclear matter 
[@, f|, ||, ||, 0. Indeed, for a wide range of model parameters such as the bag 
constant, bag models predict that the quark matter phase is absolutely stable 
i.e. its energy per baryon at zero pressure is lower than the one of 56 Fe. If this 
is true, this has important consequences in physics and astrophysics J7|] lead- 
ing e.g. to the possibility of so called "strange stars" [fj, |7j which are neutron 
stars purely consisting of quark matter in weak equilibrium with electrons. Of 
course, to check the model dependence of such findings it is important to per- 
form the corresponding calculations also in models different from the MIT bag 
model. In a recent work by Buballa and Oertel H] the equation of state (EOS) 
of quark matter was investigated in the framework of the Nambu-Jona-Lasinio 
(NJL) model with three quark flavors. Applying this model it was found that 
strange quark matter is not absolutely stable. This would rule out the exis- 
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tence of strange stars. On the other hand, the possibility of quark phases in 
the interior of neutron stars is in principle not excluded by this result - even 
though this possibility gets energetically less likely Only a detailed phase 
transition calculation can answer the question which effect the findings in || 
have on the existence of quark phases inside neutron stars. This is what we 
are aiming at in the present work. 

In principle, for the description of a neutron star which consists of a quark 
phase in its center and a surrounding hadronic phase (and, as we shall discuss 
below, a mixed phase in between) we need models for both phases. The most 
favorite case would be to have one model which can reliably describe both 
phases. So far, there are no such models. Therefore, we will use various ver- 
sions of the relativistic mean field model to parametrize the hadronic phase. 
For the quark phase we follow Buballa and Oertel || in using the three-flavor 
version of the NJL model. The NJL model has proved to be very successful in 
the description of the spontaneous breakdown of chiral symmetry exhibited by 
the true (nonperturbative) QCD vacuum. It explains very well the spectrum 
of the low lying mesons which is intimately connected with chiral symmetry as 
well as many other low energy phenomena of strong interaction [[], 10, 11 1. At 



high enough temperature and/or density the NJL model predicts a transition 
to a state where chiral symmetry becomes restored. Despite that promising 
features which at first sight might suggest the NJL model as a good candi- 
date for modeling both the low and high density region of a neutron star this 
model has one important shortcoming, namely it does not confine quarks. At 
low densities, however, the bulk properties of strongly interacting matter are 
significantly influenced by the fact that quarks are confined there. Therefore, 
we cannot expect that the NJL model gives reliable results for the EOS at low 
densities. Thus we will use the relativistic mean field model to describe the 
confined phase. At higher densities, however, the quarks are expected to be 
deconfined. There we expect the NJL model to be applicable since the lack of 
confinement inherent to this model is irrelevant in that regime. The interest- 
ing feature of the NJL model is that it reflects the chiral symmetry of QCD. 
Clearly, it would be preferable to have a Lagrangian for the hadronic phase 
which also respects chiral symmetry like e.g. the one constructed in [12| for 
the two-flavor case and the SU(3) generalizations [13|, Such Lagrangians, 
however, are more complicated to deal with. First applications to neutron 
star matter seem to indicate that the modifications are rather small as com- 
pared to the relativistic mean field models used here Hl5| . For simplicity, we 
therefore will restrict our considerations to the much simpler extensions of the 
Walecka model which include hyperonic degrees of freedom (relativistic mean 



3 



field models). 

The paper is organized as follows: In Sec. ||] we discuss how the EOS for 
the hadronic phase of a neutron star is calculated within several variants of 
relativistic mean field models. We keep brief here since such models are fre- 
quently used and well documented in the literature (cf. e.g. 0). In Sec. |3| 
we apply the NJL model to the description of the possible quark phase of the 
neutron star. Here we present much more details as compared to Sec. § since 
to the best of our knowledge it is the first time that the NJL model is applied 
to the description of the quark phase in a neutron star. Sec. [| is devoted to 
the construction of the phase transition and to the application of the complete 
EOS to the internal structure of the neutron star. Finally we summarize and 
discuss our results in Sec. ||. 



2 Hadronic matter 

Neutron stars cover a wide range of densities. From the surface of the star 
which is composed of iron with a density of e ~ 8g/cm 3 the density can in- 
crease up to several times normal nuclear matter density (eo = 140MeV/fm 3 
« 2.5 x 10 14 g/cm 3 ) in the center of the star. Since there is no single theory 
that covers this huge density range, we are forced to use different models to 
meet the requirements of the various degrees of freedom opened up at differ- 
ent densities. For subnuclear densities we apply the Baym-Pethick-Sutherland 



EOS [16]. The degrees of freedom in this EOS are nuclei, electrons and neu- 
trons. The background of neutrons appears above neutron drip density (edrip ~ 
4 x 10 11 g/cm 3 ) when the most weakly bound neutrons start to drip out of the 
nuclei which themselves get more and more neutron rich with increasing den- 
sity. For a detailed discussion of the Baym-Pethick-Sutherland EOS see also 
Hfl. We also refer to [lTj where a relativistic mean field model is extended to 



also describe this low density range. 

At densities of about normal nuclear density €q the nuclei begin to dissolve 
and merge together and nucleons become the relevant degrees of freedom in this 
phase. We want to describe this phase in the framework of the relativistic mean 
field (RMF) model which is widely used for the description of dense nuclear 



matter [[18], [T^, |2"0fl . For an introduction to the RMF model see e.g. 0. We 
use three EOS's calculated by Schaffner and Mishustin in the extended RMF 
model pg (denoted as TM1, TM2, GL85) and one by Gosh, Phatak and Sahu 



||]. For the latter one we use GPS as an abbreviation. These models include 



hyperonic degrees of freedom which typically appear at e ~ 2 — 3 eo- Table [I] 
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Table 1: Nuclear matter properties of the hadronic EOS's. The saturation 
density and the binding energy is denoted by po and B/A, the incompressibility 
by K, the effective mass by m* N /m n and the symmetry energy by a sym . The 
particle compositions are shown at the bottom of the table. 

shows the nuclear matter properties and the particle composition of the four 
EOS's. The RMF EOS's are matched to the Baym-Pethick-Sutherland EOS 
at densities of e ~ 10 14 g/cm 3 « eg. Even if the relevant degrees of freedom are 
specified (in the RMF case basically nucleons and hyperons) the high density 
range of the EOS is still not well understood. The use of different hadronic 
models should reflect this uncertainty to some degree. In the following we 
denote the phase described by the Baym-Pethick-Sutherland EOS and by the 
RMF model as the hadronic phase (HP) of the neutron star. 



3 Quark phase 



To describe the deconfined quark phase (QP) we use the Nambu-Jona-Lasinio 



(NJL) model [22] with three flavors |23[ in Hartree (mean field) approximation 
(for reviews on the NJL model cf. || |TU|, |TTJ). The Lagrangian is given by 
(cf. li) 



£ = q(i@ - m) q + G^[(q\ k q) 2 + {qi^\ k qf 



k=0 



K [ det f (q (1 + 75 ) q) + det f (q (1 - 7s ) q) 



where q denotes a quark field with three flavors, u, d, and s, and three colors. 
rh = diag(m tt , m<f , m s ) is a 3 x 3 matrix in flavor space. For simplicity we use the 
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isospin symmetric case, m u = = m q . The A& matrices act in flavor space. 
For k = 1, . . . , 8 they are the generators of SU(3) f while Ao is proportional to 
the unit matrix in flavor space (see |10] for details). The four-point interaction 
term ~ G is symmetric in SUv(3) x SUa(3) x Uy(l) x Ua(X)- In contrast, the 
determinant term ~ K which for the case of three flavors generates a six-point 
interaction breaks the {7a(1) symmetry. If the mass terms are neglected the 
overall symmetry of the Lagrangian therefore is SUv(3) x SUa{3) x t/y(l). 
In vacuum this symmetry is spontaneously broken down to SUy(3) x 
which implies the strict conservation of baryon and flavor number. The full 
chiral symmetry - which implies in addition the conservation of the axial flavor 
current - becomes restored at sufficiently high temperatures and/or densities. 
The finite mass terms introduce an additional explicit breaking of the chi- 
ral symmetry. On account of the chiral symmetry breaking mechanism the 
quarks get constituent quark masses which in vacuum are considerable larger 
than their current quark mass values. In media with very high quark densities 
constituent and current quark masses become approximately the same (con- 
cerning the strange quarks this density regime lies far beyond the point where 
chiral symmetry is restored). 

The coupling constants G and K appearing in (|l[) have dimension energy -2 
and energy -5 , respectively. To regularize divergent loop integrals we use for 
simplicity a sharp cut-off A in 3-momentum space. Thus we have at all five 
parameters, namely the current quark masses m q and m s , the coupling con- 
stants G and K, and the cut-off A. Following |2B| we use A = 602.3 MeV, 
GA 2 = 1.835, KA 5 = 12.36, m g = 5.5 MeV, and m s = 140.7 MeV. These pa- 
rameters are chosen such that the empirical values for the pion decay constant 
and the meson masses of pion, kaon and r/ can be reproduced. The mass of 
the r] meson is underestimated by about 6%. 

We treat the three-flavor NJL model in the Hartree approximation which 
amounts to solve in a selfconsistent way the following gap equations for the 
dynamically generated constituent (effective) quark masses: 

m* =nii- AG{qiqi) + 2K (qjqj) (q k q k ) (2) 

with (i, j, k) being any permutation of {u, d, s). At zero temperature but finite 
quark chemical potentials the quark condensates are given by 

Pi,<ei<A ( ' \/( m » 2 + P * I •J«) 2 + P 2 

where we have taken the number of colors to be A^ c = 3. p l F denotes the Fermi 
momentum of the respective quark flavor i. It is connected with the respective 
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quark chemical potential fi{ via 



p t F = ^-(m^e^-m*). (4) 
The corresponding quark particle number density is given by 

Pi=2iVc y (2^ = — • (5) 

Ip1<p5, 

For later use we also introduce the baryon particle number density 

P = \i.Pu + Pd + ps) ■ (6) 

The Eqs. (||), @ serve to generate constituent quark masses which decrease 
with increasing densities from their vacuum values of rn* q vac = 367.7 MeV and 
m* vac = 549.5 MeV, respectively. 

Before calculating the EOS we would like to comment briefly on the Hartree 
approximation to the NJL model which we use throughout this work. This 
treatment is identical to a leading order calculation in the inverse number 
of colors 1/N C ||10|| . In principle, one can go beyond this approximation by 
taking into account 1/N C corrections in a systematic way. This amounts in 
the inclusion of quark-antiquark states (mesons) as RPA modes in the ther- 
modynamical calculations f24| , |25| . Several things then change: First of all, 
these mesons might contribute to the EOS. We are not aware of a thorough 
discussion of such an EOS for three flavors with finite current quark masses. 
The two-flavor case is discussed in . Qualitatively the masses of the meson 
states rise above the chiral transition point. Therefore, they should become 
energetically disfavored and thus less important. An additional technical com- 
plication arises due to the fact that the relation between the Fermi energy and 
the chemical potential becomes nontrivial. Instead of one has to solve an 
additional gap equation for each flavor species. These gap equations are cou- 
pled to the gap equations for the constituent quark masses given in @. We 



refer to [10] for details. For simplicity we will restrict ourselves in the following 
to the Hartree approximation and comment on the possible limitations of that 
approach in the last section. 

Coming back to the EOS we also need the energy density and the pressure 
of the quark system. In the Hartree approximation the energy density turns 
out to be H 

Pf 

3 



. 7T 2 



€njl= E "2 dpp 2 J(m*y+p* +B eS (7) 
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while pressure and energy density are related via 



Pnjl + ^njl — E 



Pi Pi 



i=u,d,s 

where the effective bag pressure B ef[ is given by 

B e fi = Bq — B 

with 



(8) 



(9) 



B 



E 



7T- 



dpp 2 ( ^K) 2 + P 2 - \/(™i) 2 + P 2 ) ~ 2G(q iqi ) 



and 



+ 4:K(uu)(dd)(ss) 



Bn = B\ 



(10) 



,,,=, d =,.M> = (217.6 MeV) 4 . 

Note that B depends implicitly on the quark densities via the (density depen- 
dent) constituent quark masses. The appearance of the density independent 
constant Bq ensures that energy density and pressure vanish in vacuum. We 
note here that this requirement fixes the density independent part of B e g which 
influences the EOS via (|7|), (§) and therefore the possible phase transition to 
quark matter. We will come back to this point in the last section. In what 
follows we shall frequently compare the results of the three-flavor NJL model 
with the simpler MIT bag model || || . For that purpose it is important to real- 
ize that the NJL model predicts a (density dependent) bag pressure B e s while 
in the MIT bag model the bag constant -Bmit is a density independent free 
parameter. There usually also the quark masses mf lIT are treated as density 
independent quantities. (An exception is the model discussed in |^6|, |27J which 



uses density dependent effective quark masses caused by quark interactions in 
the high density regime.) In the bag model energy density and pressure of the 
quark system are given by 



p' f 



i=u,d,s 



and 



PMIT + CMIT — E Pi Pi ■ 



(12) 



(13) 



i=u,d,t 
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Figure 1: Fermi momenta of the quarks as a function of baryon particle number 
density for a charge neutral system of quarks and electrons in weak equilibrium. 
A denotes the cut-off introduced to regularize the NJL model. p = 0.17fm~ 3 
denotes nuclear saturation density. 



Suppose now that the densities are so high that in the three-flavor NJL model 
the effective quark masses have dropped down to the current quark masses. 
In this case, energy density and pressure take the form of the respective ex- 
pressions in the MIT bag model with mf llT = rrii and Burr — B . However, a 
word of caution is in order here. For very high quark particle number densities 
the corresponding Fermi momenta become larger than the momentum cut-off 
A introduced to regularize the NJL model. In this case the results of the NJL 
model become unreliable. E.g. the upper limit of the momentum integration 
in (|5|) would be no longer given by the Fermi momentum but by the cut-off 
A which would be clearly an unphysical behavior of the model. Thus for all 
practical purposes one should always ensure that in the region of interest the 
Fermi momenta are smaller than the momentum cut-off A. Fig. p] shows the 
Fermi momenta of the quarks as a function of the baryon particle number den- 
sity (for a charge neutral system of quarks and electrons in weak equilibrium; 
cf. next paragraph for details). Obviously all Fermi momenta stay below the 
cut-off A for the region of interest. We will come back to that point at the end 
of this section. 

The QP which might be found in the center of a neutron star consists of 
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u, d, and s quarks and electrons in weak equilibrium, i.e. the weak reactions 



d — > u + e~ + v e - , 

s — > u + e~ + v e - , (14) 
s + u * — > d + u 

imply relations between the four chemical potentials fj, u , fid, A*s, Pe which read 

Ps = Pd = Hu + Pe- (15) 

Since the neutrinos can diffuse out of the star their chemical potentials are 
taken to be zero. The number of chemical potentials necessary for the descrip- 
tion of the QP in weak equilibrium is therefore reduced to two independent 
ones. For convenience we choose the pair (/x n , fi e ) with the neutron chemical 
potential 

yL n = n u + 2yi, d . (16) 

In a pure QP (in contrast to quark matter in a mixed phase which we will 
discuss later) we can require the QP to be charge neutral. This gives us an 
additional constraint on the chemical potentials via the following relation for 
the particle number densities: 

2 11 

gP« - gPd - 3PS - Pe = (17) 

where p e denotes the electron particle number density. Neglecting the electron 
mass it is given by 



Pe 



f4 



3tt 2 ' 

Utilizing the relations ( [L5D and fllTD the EOS can now be parametrized by 
only one chemical potential, say /i n . At this point it should be noted that the 
arguments given here for the QP also holds for the HP. There one also ends up 
with two independent chemical potentials (e.g. fi n and /x e ) if one only requires 
weak equilibrium between the constituents of the HP and with one chemical 
potential (e.g. /x n ) if one additionally requires charge neutrality. As we will 
discuss later, the number of independent chemical potentials plays a crucial 
role in the formulation of the Gibbs condition for chemical and mechanical 
equilibrium between the HP and the QP. 

In the pure QP total energy density and pressure are given by the respective 
sums for the quark and the electron system, i.e. 

6 = 6NJL + ^ (19) 
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Figure 2: Quark condensates as a function of baryon particle number density 
for a charge neutral system of quarks and electrons in weak equilibrium. Note 
that all condensates have negative values (cf. (Q)). 

and 

V = PNJL + ^ (20) 

where the system of electrons is treated as a massless ideal gas. One obtains 
the analogous expressions for the MIT bag model if £njl and pnjl are replaced 
by the respective MIT expressions (|T2"D and (|T3|). 

Demanding weak chemical equilibrium ( |T5D and charge neutrality ( |T7D as 
discussed above all thermodynamic quantities as well as quark condensates, 
effective quark masses etc. can be calculated as a function of one chemical 
potential fi n . The curves in Fig. [1| as well as in Figs. Q-|7j which we shall discuss 
in the following are obtained by varying /i n while obeying simultaneously the 
constraints ([T^) and ([IT]). 

Figs. and § show the quark condensates and the effective quark masses, 
respectively, as a function of the baryon particle number density. Note that 
we start already at a density as high as two times nuclear saturation density 
Po = 0.17 fm -3 since we want to describe only the high density regime of the 
neutron star with quark degrees of freedom while for low densities we use the 
hadronic EOS described in the previous section. Concerning the low density 
regime of the three-flavor NJL model we refer to [|| for details. There it was 
shown that the energy per baryon of a charge neutral system of quarks and 
electrons in weak equilibrium (described by the NJL model and a free electron 
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Figure 3: Effective quark masses as a function of baryon particle number 
density for a charge neutral system of quarks and electrons in weak equilibrium. 
In addition the strange quark chemical potential /i s is plotted in the region 
where it meets the constituent strange quark mass. This marks the point 
where strange quarks appear in the system (see also Fig. [1]). 



gas) shows a minimum somewhat above two times po- This implies that in 
the density region below this minimum the pressure is negative. We are not 
interested in the (low density) part of the EOS with negative pressure since it 
cannot be realized in a neutron star. In the region of interest Figs. H and || show 
that the strange quark condensate and the effective strange quark mass stay 
constant until the strange chemical potential fi s overwhelms the strange quark 
mass. Only then according to @ the strange quark particle number density 
p s and the corresponding Fermi momentum p s F (cf. Fig. [I]) become different 
from zero causing a decrease of \(ss)\ and m*. Note that all condensates 
have negative values (cf. (|3|)). One might wonder why the dropping of the 
condensates of the light up and down quarks does not decrease the strange 
quark mass (and condensate) due to the last coupling term in (Ep. Indeed, 
strange quark mass and condensate have dropped in the low density region 
(not shown here) from their vacuum values down to the plateaus shown in 
Figs. 0, || due to their coupling to the up and down quark condensates. In the 
plateau region, however, these condensates have already decreased so much 
that their influence on the strange quark mass is diminished. We refer to || 
for details. As we shall see below, the large plateau value of the strange quark 
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Figure 4: The bag pressure (to the power of 1/4) as function of the baryon 
particle number density for the same conditions as described in Fig. || and [|. 



mass will have considerable influence on the phase structure in the interior of 
neutron stars. 

Fig. £| shows the bag pressure B ef[ as a function of the baryon particle 
number density. After staying more or less constant up to roughly 5 times 
nuclear saturation density it starts to increase towards B which, however, it 
will reach only very slowly. Again the rising of -B e fr can be traced back to the 
strange quarks which come into play at high densities. 

Thermodynamic quantities are shown in Figs. |5|-[7[ For comparison various 
curves calculated within the MIT bag model are added. The curves labeled 



with a specific value of the bag pressure are obtained from (12,13) in weak 
equilibrium where the respective value of £>mit and the current quark masses 
m q = 5.5 MeV and m s = 140.7 MeV are used. In contrast to that for the 
curve labeled with "MIT" we have used the plateau values of the bag pressure 
Emit = 170 MeV (cf. Fig. f|) and of the strange quark mass m^ IT = 465 MeV 
(cf. Fig. §). For up and down quarks we have used the current quark mass 
values also here. Fig. [5] shows the energy per baryon as a function of the 
baryon particle number density. We find that the results of the NJL model 
calculation cannot be reproduced by a bag model using the current strange 
quark mass - no matter which bag pressure is chosen. As already discussed 
above, the reason simply is that in the NJL model up to four times po there 
are no strange quarks in a system which is in weak equilibrium (cf. Fig. [1]). 
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Figure 5: Energy per baryon versus baryon particle number density for the 
NJL model and the MIT bag model for various values of the bag pressure and 
the strange quark mass. The curve labeled with "MIT" uses a bag pressure 
of -Bmit = 170 MeV and a strange quark mass of m^ IT = 465 MeV. All other 
bag model curves use the current strange quark mass of 140.7 MeV. See main 
text for details. 




On the other hand, in bag models using the much lower current strange quark 
mass one finds a reasonable amount of strange quarks already at vanishing 
pressure which typically corresponds to 2—3 times po P6fl . In contrast to that, 



a bag model with the plateau values for bag pressure and strange quark mass 
(denoted as "MIT" in the figures) yields a very good approximation to the NJL 
result for the energy per baryon up to 6—7 times Pq. For higher particle number 
densities the NJL result bends over and can be better described by bag models 
using the current strange quark mass and higher bag pressures (roughly B Q ). 
All these findings also apply to the interpretation of Fig. |6] which shows the 
total pressure of the system versus the baryon chemical potential. Comparing 
the two curves with the same bag constant labeled with "-Bmit = 170 MeV" and 
with "MIT", respectively, one observes that the latter one has a significantly 
lower pressure. This is due to the use of the much larger effective strange quark 
mass of mf IT = 465 MeV in the latter compared to the current strange 

quark mass of 140.7 MeV used in the former. The p versus fi n relation is an 
important ingredient for the construction of the phase transition from hadronic 
to quark matter inside a neutron star. We note already here, however, that 
we need in addition the thermodynamical relations also for a quark-electron 
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Figure 6: Pressure versus neutron chemical potential for the NJL model and 
the MIT bag model for various values of the bag pressure and the strange 
quark mass. See Fig. [| and main text for details. 



system away from the charge neutral configuration to describe correctly the 
phase transition (see below). 

The outlined picture concerning the comparison of NJL and bag models is 
somewhat modified when looking at Fig. [7] which shows the total pressure as 
a function of the energy density. EOS's in the form p(e) enter the Tolman- 
Oppenheimer-Volkoff [^8| equation which in turn determines the mass-radius 
relation of neutron stars. We see that in the lower part of the plotted energy 
density range the EOS in Fig. [7] is reasonably well described by MIT bag models 
with the plateau value B\[i T = 170 MeV no matter which quark masses are 
chosen (current or effective quark masses). The reason is that the p(e) relation 
is not very sensitive to the quark masses. This has already been observed in 
a somewhat different context in ]3B|]. Going to higher densities the strange 
quarks enter the game and the EOS in Fig. [7| obtained from the NJL model 
starts to deviate from the EOS of the MIT bag models with the plateau value 
-E?mit = 170 MeV. For very high densities the pressure determined from the 
NJL model becomes comparable to the one calculated in the bag model with 
a high bag constant (roughly B ). It is interesting to note that the deviation 
between the NJL curve and the "MIT" curve starts to increase in Fig. [7] much 
earlier than in Figs. [5]and|| This shows that the pressure versus energy density 
relation is much more sensitive to the detailed modeling than the relations 
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Figure 7: Pressure versus energy density for the NJL model and the MIT 
bag model for various values of the bag pressure and the strange quark mass. 
The upper scale shows the corresponding particle number densities. e = 
140 MeV/fm 3 . See Fig. [5] and main text for details. 



shown in Figs. |5] and |. 

Before constructing the phase transition inside the neutron star let us 
briefly discuss the limitations of the NJL model in the form as we have treated 
it here. As a typical low energy theory the NJL model is not renormaliz- 
able. This is not an obstacle since such theories by construction should be 
only applied to low energy problems. In practice the results depend on the 
chosen cut-off or, to turn the argument around, the NJL model is only prop- 
erly defined once a cut-off has been chosen. This cut-off serves as a limit for 
the range of applicability of the model. Here we have used one cut-off A for 
the three-momenta of all quark species. Concerning the discussion of other 



cut-off schemes and their interrelations we refer to [0, [TTJ] . When the density 
in the quark phase gets higher the Fermi momenta of the quarks rise due to 
the Pauli principle. Eventually they might overwhelm the cut-off of the NJL 
model. At least beyond that point the model is no longer applicable. We have 
made sure in our calculations that this point is never reached (cf. Fig. [I]). In 
addition, at very high densities one presumably enters a regime which might 
be better described by (resummed) perturbation theory While the nonpertur- 
bative features of the NJL model vanish with rising density, medium effects 
as mediated e.g. by one-gluon exchange grow with the density [^6], |27f . To 
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summarize, concerning the calculation of the EOS it turns out that the NJL 
model should neither be used at low densities where confinement properties 
are important nor at very high densities where the NJL model as a low energy 
theory leaves its range of applicability. However, the NJL model might yield 
reasonable results in a window of the density range where confinement is no 
longer crucial but chiral symmetry as a symmetry of full QCD remains to be 
important. 

4 Phase transition and neutron stars 

In the previous sections we have discussed the underlying EOS's thought to 
reflect the properties of confined hadronic matter (HP) and deconfined quark 
matter (QP) in its particular regime of applicability. Applying these EOS's 
we want to calculate in this section the phase transition from the HP to the 
QP to see which phase is the favored one at which densities. (The existence 
of a QP inside the neutron star of course requires the phase transition density 
to be smaller than the central density of the star.) 

It is worth to point out which phase structure is in principle possible if a 
hadronic model and the NJL model are connected at a certain density value 
Pdeconf (which is dynamically determined in the present work by a Gibbs con- 
struction as we shall discuss below). At density pdeconf we assume a first order 
phase transition from confined hadronic to deconfined quark matter. Even 
without a matching to a hadronic model the NJL model already exhibits a 
transition, namely from a low density system with broken chiral symmetry 
to a high density system where chiral symmetry is restored. The respective 
density is denoted by p c hirai- For densities larger than p c hirai the Goldstone 
bosons which characterize the chirally broken phase are no longer stable but 
can decay into quark- ant iquark pairs. If p c hirai was larger than pdeconf the fol- 
lowing scenario would be conceivable: There would be three phases, namely (i) 
a hadronic, i.e. confined phase at low densities, (ii) a phase where quarks are 
deconfined but massive (in this phase e.g. pions would still appear as bound 
states), and (iii) a high density phase where quarks are deconfined and their 
masses are so low that all mesons can decay into quarks. Had we neglected all 
current quark masses, the quarks in the third phase would be exactly massless. 
With finite current quark masses, however, the constituent quark masses keep 
on dropping with rising density in the third phase (cf. Fig. ^) . This definitely 
interesting scenario with three phases is not realized in our model. It turns 
out that the deconfinement phase transition happens far beyond the chiral 
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transition, i.e. p c hirai < Pdeconf- Therefore, only the phases (i) and (iii) appear 
here. 

In principle, since we assume the deconfmement phase transition to be of 
first order these two phases can coexist in a mixed phase. Indeed, it was first 
pointed out by Glendenning that beside a HP and a QP also this mixed phase 
(MP) of quark and hadronic matter may exist inside neutron stars [0, |29|. 
(For a discussion of the geometrical structure of the MP and its consequences 
for the properties of neutron stars see 0.) This possibility was not realized in 
previous calculations due to an inadequate treatment of neutron star matter as 
a one-component system (one which can be parameterized by only one chemical 
potential). As we have already discussed, the treatment of neutron star matter 
as a charge neutral phase in weak equilibrium indeed reduces the number of 
independent chemical potentials to one. But the essential point is that - if a 
MP exists - charge neutrality can be achieved in this phase e.g. with a positively 
charged amount of hadronic matter and a negatively charged amount of quark 
matter. Therefore it is not justified to require charge neutrality in both phases 
separately. In doing so we would "freeze out" a degree of freedom which in 
principle could be exploited in the MP by rearranging electric charge between 
both phases to reach "global" charge neutrality. A correct treatment of the 
phase transition therefore only requires both phases to be in weak equilibrium, 
i.e. both phases still depend on two independent chemical potentials. We have 
chosen the pair (/i n , /i e ). Such a system is called a two-component system. The 
Gibbs condition for mechanical and chemical equilibrium at zero temperature 
between both phases of the two-component system reads 

PHp{Hn, AO = PQp{Hn, He) = PMP- (21) 

Using Eq. fl2T|) we can calculate the equilibrium chemical potentials of the 
MP where pup = Vqp holds. Fig. § illustrates this calculation. The HP— >MP 
phase transition takes place if the pressure of the charge neutral HP (white line) 
meets the pressure surface of the QP (NJL). Up to this point the pressure of 
the QP is below the pressure of the HP making the HP the physically realized 
one. At higher pressure the physically realized phase follows the MP curve 
which is given by the Gibbs condition (EH]). Finally the MP curve meets the 
charge neutral QP curve (white line) and the pressure of the QP is above the 
pressure of the HP, making the QP the physically realized one. For every point 
on the MP curve one now can calculate the volume proportion 



Vqp + Vhp 
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Figure 8: Gibbs phase construction of a two-component system. Plotted is the 
pressure surface of the hadronic phase (RMF) and of the quark phase (NJL) as 
a function of the two independent chemical potentials /i n , /jL e . EOS of the HP 
is GPS, EOS of the QP uses the NJL model. The white lines HP and QP show 
the pressure of the hadronic and the quark phase under the condition of charge 
neutrality. At low pressure matter is in its charge neutral HP. The intersection 
curve MP corresponds to the mixed phase. This curve is the solution of the 
Gibbs condition ([H]). At very high pressure matter consists of a pure QP. 
Also shown are the central pressures of a typical M = 1.4 M & neutron star 
and of an M = 1.6 M© star. (M denotes the mass of the sun.) Clearly 
neither achieves a central pressure necessary to undergo a phase transition to 
deconfined matter. 
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occupied by quark matter in the MP by imposing the condition of global charge 
neutrality of the MP 

XP? P + (l-x)pf P = 0. (23) 

Here pf p and p^ p denote the respective charge densities. From this, the energy 
density e of the MP can be calculated by 

£mp = X egp + (1 - X) e np- (24) 

Along the MP curve the volume proportion occupied by quark matter is mono- 
tonically increasing from x = to x = 1 where the transition to the pure QP 
takes place. 

Taking (i) the charge neutral EOS of the HP at low densities (Sec. |2|), (ii) 
Eq. (0), (||), and (||) for the MP, and (iii) the charge neutral EOS of the 
QP (Sec. we can construct the full EOS in the form p = p(e). For simplicity 
we denote this EOS as the hybrid star EOS. Fig. [| shows this EOS if we apply 
GPS for the HP EOS. The hybrid star EOS consists of three distinct parts. 




5 10 15 20 

e/fo 



Figure 9: EOS in the form pressure versus energy density using GPS for the 
hadronic EOS. The shaded region (MP) corresponds to the mixed phase part of 
the EOS. eo = 140MeV/fm 3 . The mass scale in the plot shows the respective 
neutron star mass if the corresponding energy density is regarded as the central 
density of the neutron star. Obviously around the typical neutron star mass 
of 1.4 M & the neutron stars consist solely of hadronic matter. M & denotes the 
mass of the sun. 

At low densities (e^7e ) matter is still in its confined HP. At (e ~ 7e ) the 
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first droplets of deconfined quark matter appear. Above this density matter is 
composed of a mixed phase of hadronic and quark matter. This MP part of 
the EOS is shaded gray. Only at unaccessible high densities (e^l8eo) matter 
consists of a pure QP. The preceding statements refer to the use of GPS for the 
EOS of the HP. Concerning all the other variants of RMF used here (TM1, 
TM2, GL85) we have found that the HP^MP transition does not appear 
below e ~ 10 6q. As we will discuss below such high energy densities cannot be 
reached inside a stable neutron star which is described by one of these EOS. 

At this point we should note the essential difference between the treatment 
of neutron star matter as a one- and a two-component system (cf. f29|). While 
the former one leads to the well known phase transition with a constant pres- 
sure MP (like in the familiar liquid-gas phase transition of water) , we can see in 
Fig. |^ that the pressure is monotonically increasing even in the MP if we apply 
the correct two-component treatment. This has an important consequence on 
the structure of the neutron star. Since we know from the equations of hydro- 
static equilibrium - the Tolman-Oppenheimer-Volkoff (TOV) equations [E3[] 
that the pressure has to increase if we go deeper into the star, a constant pres- 
sure MP is strictly excluded from the star while a MP with increasing pressure 
can (in principle) occupy a finite range inside the star. 

To see if the densities inside a neutron star are high enough to establish a 
MP or a QP in its center we have to solve the TOV equations with a specified 
hybrid star EOS following from our phase transition calculation. From the 
solutions of the TOV equations we get a relation between the central energy 
density (or central pressure) and the mass of the neutron star (cf. Fig. |i~0|). 
The maximum possible central energy density (the critical energy density e cr i t ) 
is reached at the maximum mass that is supported by the Fermi pressure of 
the particular hybrid star EOS. Above this critical density the neutron star 
gets instable with respect to radial modes of oscillations |2j . We have applied 
the four HP EOS's (denoted by GPS, TM1, TM2 and GL85) to calculate the 
four corresponding hybrid star EOS's. (The one for GPS is shown in Fig. [].) 
We found that in no EOS the central energy density of a typical M « 1.4 M & 
neutron star is large enough for a deconfinement phase transition. (Here M Q 
denotes the mass of the sun.) The corresponding neutron stars are purely 
made of hadronic matter (HP). In Fig. § where GPS is used also the neutron 
star masses are shown as a function of the central energy density. There the 
central energy density of a M = 1.4 M & is about e ~ 3 eo which is clearly below 
e « 7 €q which is at least necessary to yield a MP core. (This is also shown in 
the context of the Gibbs construction in Fig. |8] where the central pressure of a 
M = 1.4 M Q and of a M = 1.6 M© neutron star is marked.) In Fig. |] we can 
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Figure 10: Mass of the neutron star as a function of the central energy density. 
Above a central energy density of e cr i t ~ 9 e where the maximum mass of 
the neutron star is reached the system becomes instable with respect to radial 
modes of oscillations. The shaded region (MP) corresponds to stars which 
possess a mixed phase in their center. The hadronic part of the EOS uses 
GPS. e = 140MeV/fm 3 . M denotes the mass of the sun. 

see that only near the maximum mass of M m 1.72 M & neutron stars with a 
MP core are possible. This, however, only holds for the GPS hybrid star EOS 
and only in a quite narrow mass range from M ~ 1.7—1.72 M . In the density 
range up to the critical density all other EOS's (TM1, TM2, GL85) do not show 
a phase transition at all. The critical energy densities for these EOS's are in 
the range of e crit ~ 5 — 6e while the densities for the HP^MP transition are 
above e ps 10 e . (The corresponding maximum masses are M ps 1.5— 1.8 M .) 
From this we conclude that within the model constructed here the appearance 
of deconfined quark matter in the center of neutron stars turns out to be very 
unlikely. 

5 Summary and conclusions 

We have studied the possible phase transition inside neutron stars from con- 
fined to deconfined matter. For the description of the quark phase we have 
utilized the NJL model which respects chiral symmetry and yields dynamically 
generated quark masses via the effect of spontaneous chiral symmetry break- 
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ing. We found that the appearance of deconfined quark matter in the center 
of a neutron star appears to be very unlikely, for most of the studied hadronic 
EOS's even impossible. The ultimate reason for that effect is the high value 
of the effective strange quark mass which turns out to be much higher than 
its current mass value in the whole relevant density range (cf. Fig. |3|). This 
finding, of course, is based on several assumptions which need not necessarily 
be correct. In lack of an EOS based on a full QCD calculation at zero tem- 
perature and finite nuclear density, we had to rely on simpler models for the 
EOS in different density regimes. 

Concerning the low density regime we have used various relativistic mean 
field (RMF) models. To some degree the use of different variants of the RMF 
model should reflect the uncertainties of this approach. These models are 
generalizations of the Walecka model [ 30 1 which describes the hadronic ground 
state of nuclear matter at density p quite successfully. At somewhat higher 
densities the used RMF models deal with hyperons as additional degrees of 
freedom. Since these RMF models do not have any explicit quark degrees 
of freedom we expect them to become unreliable at high densities where the 
confinement forces are screened and the hadrons dissolve into quarks. 

To describe this high density regime we have utilized the Nambu-Jona- 
Lasinio (N JL) model in its three-flavor extension. The merits of the NJL model 
are (at least) twofold: For the vacuum case, it gives a reasonable description 
of spontaneous chiral symmetry breaking and of the spectrum of the low lying 
mesons. For sufficiently high density and/or temperature, the NJL model 
exhibits the restoration of chiral symmetry. A shortcoming of the NJL model 
is that it does not confine quarks, i.e. there is no mechanism which prevents 
the propagation of a single quark in vacuum. Therefore in an NJL model 
calculation the quarks significantly contribute to the EOS also at low densities.^ 
This was the ultimate reason why we considered the NJL model only in the 
high density regime where confinement is supposed to be absent anyway while 
utilizing the hadronic RMF models to describe the confined phase. On the 
other hand, we should recall that energy density and pressure of the NJL model 
were determined such that both vanish at zero density, i.e. in a regime where 
we have not utilized the NJL model afterwards. This procedure determines 
the effective bag pressure B e g given in @ by fixing B (|TTD to (217.6 Me V) 4 . 
Clearly, this procedure is somewhat unsatisfying since the effective bag pressure 
-B e ff influences the EOS and therefore the onset of the phase transition. Indeed, 



1 Actually in the mean field approximation the quarks are the only degrees of freedom 
which contribute. For a generalization to include meson states as RPA modes in the NJL 
EOS see \M EHl. 
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if we reduce Bq by only 5—10% by hand from its original value of (217.6 MeV) 4 
we already observe drastic changes in the phase structure of the neutron star 
favoring deconfmed quark matter. On the other hand, the physical requirement 
that any model should yield vanishing energy density and pressure in vacuum 
is the only way to uniquely determine the EOS of the NJL model without any 
further assumptions. 

Possible alternatives to the use of the NJL model for the description of the 
deconfmed quark matter are the MIT bag model and the extended effective 
mass bag model |^6[ ^7fl. The latter includes medium effects due to one-gluon 
exchange which rise with density (while the effective masses of the NJL model 
decrease). As already discussed at the end of Sec. |3] for the regime of very 
high densities such a resummed perturbation theory might be more adequate. 
Furthermore, MIT bag models can be very useful in interpreting more involved 
models like the NJL model in terms of simple physical quantities like the bag 
constant and the quark masses. Therefore we have frequently compared our 
NJL model results with the MIT bag model in Sec. |3|. For a further discussion 
of the MIT bag model and its application to neutron stars see 

The distinct feature of the NJL model is that nonperturbative effects are 
still present beyond the phase transition point. It is reasonable to consider 
such effects since it was found in lattice calculations |32], [33[ that for QCD at 
finite temperature the EOS beyond the phase transition point can neither be 
properly described by a free gas of quarks and gluons nor by QCD perturba- 
tion theory |35|. Presumably this holds also for the finite density regime. 
On the mean field level the most prominent nonperturbative feature of the 
NJL model which remains present beyond the phase transition point is the 
constituent strange quark mass which is much larger than the current strange 
quark mass in the whole relevant density regime (cf. Fig. [3]). This high strange 
quark mass has turned out to be crucial for the phase transition. As we have 
shown above the EOS of the QP can be reasonably well approximated by an 
MIT bag EOS up to 5 times p using a comparatively low bag constant of 
I?mit = 170 MeV and an effective strange quark mass of m^ IT = 465 MeV. 
These are the plateau values of the corresponding quantities in the NJL model 
calculations shown in Figs. |3| and |j. In Figs. [5|-[7| the curves labeled by MIT 
use these values for the bag constant and the effective strange quark mass. 
Using such a bag constant in connection with the current strange quark mass 
would allow the existence of a QP inside a neutron star @, |3l|. This, however, 
does not remain true once a much higher effective strange quark mass is used. 
Qualitatively, the chain of arguments is that a higher mass leads to a lower 
pressure (cf. Fig. 0). This disfavors the quark phase in the Gibbs construction, 
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i.e. shifts the phase transition point to higher densities. This is the reason 
why in our calculations the existence of quark matter in the center of a neu- 
tron star is (nearly) excluded. Especially for typical neutron stars with masses 
M ~ 1.4M Q the central energy density is far below the deconfinement phase 
transition density (cf. Fig. |7|). This finding is independent of the choice of the 
version of the RMF model. This suggests that it is the NJL model with its 
large strange quark mass which defers the onset of the deconfinement phase 
transition rather than the modeling of the hadronic phase. Throughout this 
work we have used the NJL parameter set of |23|| . We have also explored the 



set given in f36fl with A = 631 MeV, GA 2 = 1.830, Kk b = 9.19, m q = 5.5 MeV, 



and m s = 138 MeV. The results are very similar to the ones presented here. 

For simplicity we have treated in the present work the NJL model in the 
Hartree approximation. In principle, going beyond the mean field approxima- 
tion might influence the order of the chiral phase transition (for related work 
towards that direction for the two-flavor case cf. [[37J). If it turned out that this 
would result in a strong first order phase transition then the effective strange 
quark mass might change more drastically and in the region of interest would 
be perhaps much lower than in the case studied in the present work. This 
would favor the appearance of quark matter in the interior of neutron stars. 
Clearly, it would be interesting to study how a more involved treatment of the 
NJL model beyond the Hartree approximation would influence our findings 
presented here. This, however, is beyond the scope of the present work. 
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